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Abstract 

In 2001, the author proposed the ultra second quantization method. The ultra 
second quantization of the Schrodinger equation, as well as its ordinary second 
quantization, is a representation of the iV-particle Schrodinger equation, and this 
means that basically the ultra second quantization of the equation is the same as 
the original iV-particle equation: they coincide in 3iV-dimensional space. 

We consider a short action pairwise potential V{xi — Xj). This means that as 
the number of particles tends to infinity, N — > do, interaction is possible for only a 
finite number of particles. Therefore, the potential depends on N in the following 
way: V/v = V((xi — x^N 1 ' 3 ). If V(y) is finite with support f2y, then as N — > oo 
the support engulfs a finite number of particles, and this number does not depend 
on N. 

As a result, it turns out that the superfluidity occurs for velocities less than 
min(A cr it, ^r); where A cr i t is the critical Landau velocity and R is the radius of the 
nanotube. 

1. In order to distinguish the classical theory in its modern understanding from the 
quantum theory, it is necessary to modify (somewhat) the ideology habitual to physicists, 
for whom the classical theory is simply the whole body of physics as it existed in the 19th 
century before the appearance of quantum theory. Actually, the correct meaning is that 
the classical theory is the limit of the quantum one as h — > 0. 

Thus, Feynman correctly understood that spin is a notion of classical mechanics. In- 
deed, it is obtained via a rigorous passage from quantum mechanics to classical mechanics 
[1]. In a similar same way, the polarization of light does not disappear when the frequency 
is increased, and is therefore a property of geometric rather than wave optics, contrary to 
the generally accepted belief, which arose because the polarization of light was discovered 
as the result of the appearance of wave optics. 

Consider a "Lifshits hole", i.e. a one-dimensional Schrodinger equation with potential 
symmetric with respect to the origin of coordinates with two troughs. Its eigenfunctions 
are symmetric or antisymmetric with respect to the origin. As h — > this symmetry 
remains, and since the square of the modulus of the eigenfunction corresponds to the 
probability of the particle to remain in the troughs, it follows that in the limit as h — > 0, 
i.e., in the "classical theory", for energies less than those required to pass over the barrier, 
the particle is simultaneously located in two troughs, although a classical particle cannot 
pass through the barrier. Nevertheless, this simple example shows how the ideology of 
the "classical theory" must be modified. 

To understand this paradox, one must take into consideration the fact that the sym- 
metry must be very precise, up to "atomic precision" , and that stationary state means a 
state that arises in the limit for "infinitely long" time. 

When we deal with nanotubes whose width is characterized by "atomic" or "quantum" 
dimensions, then new unexpected effects occur in the "classical" theory. Thus, already 
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in 1958 [2], I discovered a strange effect of the standing longitudinal wave type in a 
slightly bent infinite narrow tube, for the case in which its radius is the same everywhere 
with atomic precision. It was was impossible at the time to implement this effect in 
practice, which would have allowed to obtain a unimode laser, despite A.M.Prokhorov's 
great interest in the effect. 

2. Now let us discuss the notion known as "collective oscillations" in classical physics 
and as "quasiparticles" in quantum physics. In classical physics, it is described by the 
Vlasov equation for selfcompatible (or mean) fields, in quantum physics, by the Hartrey 
(or the Hartrey-Fock) equation. 

(1) Variational equations depend on where (i.e., near what solutions of the original 
equation) we consider the variations. For example, in [3, 4, 5] we considered variations 
near a microcanonical distribution in an ergodic construction, while in [7, 8, 9, 10] this was 
done near a nanocanonical distribution concentrated on an invariant manifold of lesser 
dimension, i.e., not on a manifold of constant energy but, for example, on a Lagrangian 
manifold of dimension coinciding with that of the configuration space. 

(2) Let us note the following crucial circumstance. The solution of the variational 
equation for the Vlasov equation does not coincide with the classical limit for variational 
equations for the mean field equations in quantum theory. 

Consider the mean field equation in the form 

ih^<p%x)=(-^A + W t (x))<p t (x), W t (x) = U(x) + J V(x,y)\ V \y)\ 2 dy, (1) 

with the initial condition (p\t=o = ¥>o, where satisfies <po € W^°(R U ), f dx\(p (x)\ 2 = 1. 

In order to obtain asymptotics of the complex germ type [11] one must write out the 
system consisting of the Hartrey equation and its dual, then consider the corresponding 
variational equation, and, finally, replace the variations Sip and Sip* by the independent 
functions F and G. For the functions F and G, we obtain the following system of equa- 
tions: 

.dF\x) f ( 5 2 H , 5 2 H \ 

.dG'(x) t I PH , 6 2 H r , u \ 

-'— = J dy {s^mj F (v) + m^vw g ( 4 

The classical equations are obtained from the quantum ones, roughly speaking, by 
means of a substitution of the form ip = x e ^ S (the VKB method), ip* = x* e ^ S * > where S = 
S*, x = x(x,t)eC°°, S = S(x,t)eC°°. 

To obtain the variational equations, it is natural to take the variation not only of 
the limit equation for x and X* \ but also for the functions S and S*. This yields a new 
important term of the equation for collective oscillations. 

Let us describe this fact for the simplest example, which was studied in N.N.Bogolyubov's 
famous paper concerning "weakly ideal Bose gas" [12]. 

Suppose U — in equation (1) in a three-dimensional cubical box of edge L, the 
wave functions satisfying the periodicity condition (i.e., the problem being defined on the 
3-torus with generators of lengths L, L, L). Then the function 

<p(x) = L- 3/2 e i/h{px - m) , (3) 
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where p = 2nn/L, n - is an integer, satisfies the equation (1) for 
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tt = ^- + L- 3 I dxV(x). (4) 

Consider the functions F^ x \x) and G^(x), where A = 27rn/L,n ^ 0, n is an integer, 
given by 



G^\x) = L- 3 / 2 a x eW- p+x > + ^ +n W; (5) 



here 



2m 2m 

2 l„ |2_i t7 _ r-3 / j„TA/„\„iAx 



|a A | 2 -|p A | 2 = l, Va = L- j / cfa^e* 
From the system (6), we find 




p x = - P x + d — + y A -y A 2 . (7) 



In this example w = e^' 1 ' 4 ', -u* = e ^ , where s(x, t, ) = px + /3t, while the variation 
of the action for the vector (5u, 5u*^j equals \x ± Qt. 

Under a more accurate passage to the limit, we obtain V\ — > Vo = L~ 3 f dxV(x). 

Thus, in the classical limit, we have obtained the famous Bogolyubov relation (7). 
In the case under consideration u(x) = and, as in the linear Schrodinger equation, the 
exact solution coincides with the quasiclassical one. In the paper [10], the case u(x) ^ is 
investigated, and it turns out that the relation similar to (7) is the classical limit as h — > 
of the variational equation in this general case. The curve showing the dependence of f3\ 
on A is known as the Landau curve and determines the superffuid state. The value A cr 
for which superfluidity disappears is called the Landau critical level. Bogolyubov explains 
the superfluidity phenomenon in the following terms: "the 'degenerate condensate' can 
move without friction relatively to elementary perturbations with any sufficiently small 
velocity" [10], p. 210. 

However, there is no Bose-Einstein condensate whatever in these mathematical con- 
siderations, it is just that the spectrum defined for A < A cr is a positive spectrum of 
quasiparticles. This means it is metastable (see [13]). The Bose-Einstein condensate is 
not involved here, it is only needed only to show that it would be wrong to believe that 
this argument works for a classical liquid, as one might think from the considerations 
above. 

Indeed, for example, the molecules of a classical nondischarged liquid rule, 
Bose particles. For such a liquid, one can write out the iV-particle equation, having in 
mind that each particle (molecule) is neutral and consists of an even number / of neutrons. 
Thus each ith particle is a point in 3(2k + Z)-dimensional space, where k is the number of 
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electrons, Xi e R 6k+Sl , depends on the potential u(x-i), Xi G R ek+31 and we can consider 
the iV-particle equation for Xi,i = 1, . . . , N, with pairwise interaction V(xi — Xj). 

3. However, there is a purely mathematical explanation of this paradox. The thing 
is that Bogolyubov found only one series of of ppoints in the spectrum of the many par- 
ticle problem. Landau wrote "N.N. Bogolyubov recently succeeded, by means of a clever 
application of second quantization, in finding the general form of the energy spectrum of 
a Bose-Einstein gas with weak interaction between the particles" ([14], p. 43). But this 
series is not unique, i.e., the entire energy spectrum was not obtained. 

In 2001, the author proposed the method of ultra second quantization [15]; see also 
[16], [17], [18], [19], [20], 

The ultra second quantization of the Schrodinger equation, as well as its ordinary 
second quantization, is a representation of the iV-particle Schrodinger equation, and this 
means that basically the ultra second quantization of the equation is the same as the 
original iV-particle equation: they coincide in 3iV-dimensional space. However, the re- 
placement of the creation and annihilation operators by c-numbers, in contrast with the 
case of second quantization, does not yield the correct asymptotics, but it turns out that it 
coincides with the result of applying the Schroeder variational principle or the Bogolyubov 
variational method. 

For the exotic Bardin potential, the correct asymptotic solution coincides with the one 
obtained by applying the ultra second quantization method described above. In the case 
of general potentials, in particular for pairwise interaction potentials, the answer is not the 
same. Specifically, the ultra second quantization method gives other asymptotic series of 
eigenvalues corresponding to the iV-particle Schrodinger equation, and these eigenvalues, 
unlike the Bogolyubov ones (7), are not metastable. 

It turns out that the main point is not related to the Bose-Einstein condensate, but 
has to do with the width of the capillary (the nanotube) through which the liquid flows. 
If we consider a liquid in a capillary or a nanotube of sufficiently small radius the velocity 
corresponding to metastable states is not small. Hence at smaller velocities the flow will 
be without friction. 

The condition that the liquid does not flow through the boundary of the nanotube is 
a Dirichlet condition. It yields a standing wave, which can be regarded as a pair particle- 
antiparticle: a particle with momentum p orthogonal to the boundary of the tube, and 
an antiparticle with momentum —p. 

We consider a short action pairwise potential V(xi — Xj). This means that as the 
number of particles tends to infinity, N — > oo, interaction is possible for only a fi- 
nite number of particles. Therefore, the potential depends on N in the following way: 
Vn = V((xi — Xj)N 1 ^). If V(y) is finite with support Qy, then as iV — > oo the support 
engulfs a finite number of particles, and this number does not depend on N. 

As the result, it turns out that for velocities less than min(A cr , ^r)' where A cr is the 
critical Landau velocity and R is the radius of the nanotube, superfluidity occurs. 

Now let me present my own considerations, which are not related to the mathematical 
exposition. Viscosity is due to the collision of particles: the higher the temperature, 
the greater the number of collisions. In a nanotube, there are few collisions, and only 
with the walls, and those are taken into account by the author's series. It is precisely 
this circumstance, and not the Bose-Einstein condensate, which leads to the weakening 
of viscosity and so to superfluidity. What I am saying is that the main factor in the 
superfluidity phenomenon, even for liquid helium 4, is not the condensate, but the presence 
of an extremely thin capillary [21], [22]. 
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1 Ultrasecondary quantization 



In the papers [15, 23, 24, 25] the notion of ultrasecondary quantization was introduced 
for problems of quantum mechanics and statistical physics. Let us recall the notation 
and the main facts in the case of the quantization of pairs (particle-particle numbers) 
and for pairs consisting of two particles. Quantization of pairs will allow us to take into 
account pairwise correlations of particles in the construction of the asymptotics. The space 
of ultrasecondary quantization is J 7 , the Fock boson space [26], b + (x, s) is the creation 
operator of particles with number s, b~(x, s) is the annihilation operator of particles with 
number s in the space T [26], B + (x, x') is the creation operator of pairs particles, B~(x, x') 
is the annihilation operator of pairs of particles in this space. These operators satisfy the 
following commutation relations: 



[b-(x, s),b + (x>, s')] = 5 ss ,5(x - x'), \b±(x, s),b ± (x', s')] = 0, 
[B~ (xi, x 2 )B + (x'n x' 2 )] = 5(xi - x' 1 )5(x 2 - x 2 ), 

[B ± (x 1 ,x 2 ),B ± (x' 1 ,x' 2 )}=0, (8) 
[b ± (x,s),B ± (x' 1 ,x' 2 )] = [b ± (x,s),P(x' 1 ,x' 2 )]=0. 



Further, $ is the vacuum vector in the space J 7 , which possesses the following prop- 
erties: 



r(x,s)$ = 0, B-(x 1 ,x 2 )% = 0. (9) 

The variable x lies on the 3-dimensional torus of size Lx Lx L, which we shall denote by 
T. The variable s is discrete, s = 0, 1, . . . ; s is called the number or the statistical spin. 
Any vector $ of the space T can be uniquely represented in the form: 



$ = Y m Y ■■ • Y / • • • / dxi ■■■ dXkd y^ ■ ■ ■ d y™ x 

k=0 M=0 ' ' si=0 s k =0 J J 

x$ fc ,M(^i,si; . . .;a:fc,s fc ;yi,j/2; • • • ;j/2M-i,J/2jif) x (10) 
xb + (x u b + (x k ,s k )B + (y 1 ,y 2 ) B + (y 2M _ 1: y 2 M)$o, , 

where the function $k,M(%i, s i] ■ ■ ■ ] x k, s k] Vi,V2', ■ ■ ■] Vim-x, Vim) is symmetric with respect 
to transpositions of pairs of variables (xj,Sj) and (xj,Sj) and symmetric with respect to 
transpositions of pairs of variables (y2j-i,V2j) an d (2/2i-i, 2/2i)- 

In the boson case, we introduce the space J 7 ^^™ consisting of the vectors $ for 
which $fc',M' = whenever (k',M') ^ (k,M), while §k,M is a symmetric function of 
the variables x±, x 2 , ■ ■ ■ , x k , yi, y2, ■ ■ ■ , 2/2M- In the fermion case we similarly introduce the 
space ^f^ mm consisting of the vectors $ such that $k',M' — for (k',M') ^ (k,M), 
and $k,M is an antisymmetric function of the variables xi, x 2 , ■ ■ ■ , x^, yi, y2, ■ ■ ■ , 2/2M- The 
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orthogonal projector of the space T on the subspace * s °f the f° rm [15]- [25]: 

^ OO CO p « 

n?M m = Y " ' Y ■■ dx i--- dx kdyi ■ ■ ■ dy 2M x 

x6 + (xi,si) b + (x k , s k )B + (y 1} y 2 ) B + (y 2M -i,y 2 M) x 

x Symm (& _ (xi,si) b~(x k ,s k )B~(y 1 ,y 2 ) ■ . . . 

xi...x k yi...y 2 M ^ 

•B~(y 2M -i,y 2 M)^ x (11) 
x exp ( ~ / dx b + ( x i s)b~(x, s) - J J dydy'B + (y, y')B~(y, y'^j , , 

where Symm a . l _ x yi _j, 2M is the symmetrization operator in the variables Xi,...,x k , 

V\i ■ ■ ■ iViMi while the operators b + (x, s), b~(x,s), B + (y,y'), B~(y,y') are ordered in 
the Vick way [26]. The orthogonal projector of the space T on the subspace ^f^ mm is 
of the form [15]: 

gAsymm = __ £\ . . / . . . J dx ± . . . dx k d Vl . . . dy 2M X 

si=0 s k =0 

xb + (x 1 , si) & + (^, s k )B + (y 1 , y 2 ) B + (y 2M _ 1 , y 2M ) x 

x Asymm (6~(xi,Si) b~(x k , s k )B~(y 1 , y 2 ) ■ . . . 

x 1 ...x k y 1 ...y 2M V 



■B~(y 2M -i,y2M) S ) x 

x exp ( ~ Y J dx b + ( x , s)b~(x, s) - J J dydy'B + (y, y')B~(y, y')^j , 

where Asymm a . 1 __ M is the antisymmetrization operator in the variables x±, . . . ,x k , 

y\i ■ ■ ■ iyiM- Here and everywhere else below, unless explicitely stated, the operators 
b + (x, s), b~(x, s), B + (y,y'), B~(y,y') are ordered in the Vick way. 

Further, we consider the system of N identical particles on the torus T. We shall 
assume that the Hamiltonian for N bosons or N fermions has the form: 

t2 N N N 

^ = -^E A ^EE%-^ ( 12 ) 

3=1 j=l l=j+l 

According to [15], in the boson case, to this operator corresponds the ultra secondary 
quantized Hamiltonian 



h b = Y T[m\ Y ■ ■ ■ Y / • • • / dxi ■■■ dXkd y i ■ ■ ■ d y™ x 

fc=0 M=0 ' ' si=0 s k =0 J J 

xb + (x 1 , si) b + (x k , s k )B + (y 1} y 2 ) B + (y 2M _ u y 2M )H k+2M x 

x Symm (&~(a;i,si) b~(x k ,s k )B~(y 1 ,y 2 ) ■ . . . 

x 1 ...x k y 1 ...y 2 M ^ 

■B~(y 2M -i,y2M)^ x (13) 
x exp ( ~~ Y J dx b + ( x , s)b~(x, s) - J J dydy'B + (y, y')B~(y, y'^j , , 



while in the fermion case the corresponding operator is Hp and is expressed by a similar 
formula in which Symm is replaced by Asymm. As in (12) and (13), the ultrasecondary 

quantized operator A is assigned [15] to any iV-particle operator 



l l 

2 2.8 d 



A N \ xi,..., x N ; -i - — , 



dxi Oxn 

For example, to the unit operator in the boson case, we assign the following ultrasec- 
ondary quantized unit operator: 



e b = J2 Y ■ ■ ■ Y / • • • / dxi ■ ■ ■ dXkd y i ■ ■ ■ d y™ x 

fc=0 M=0 ' ' si=0 s k =0 J J 

xb + (xi,s l ) b + (x k ,s k )B + (y 1 ,y 2 ) B + (y 2M -i,y2M) x 

x Symm (6~(a;i,si) b~(x k ,s k )B~(y 1 ,y 2 ) ■ . . . 

x 1 ...x k y 1 ...y 2M ^ 

•B~(y 2M -i,y2M)j x 
x exp ( ~ Y J dx b + ( x , s)b~(x, s) - J J dydy'B + (y, y')B'{y, y')^j , , 

which is the sum of the projectors (11). Similarly, in the fermion case, the ultrasecondary 
quantized unit operator is 



oo oo 

r Asymm 
k,M > 

fc=0 M=0 



with Symm replaced by Asymm in formula (14). 
Consider the following eigenvalue problem 

H B>F $ = \E B $, E<$> ^ 0, (14) 
in the Bose and Fermi cases. The following assertion, proved in [15] is valid: in the 
subspaces F^m™ an d °f ^ e s V ace the operators Hb and Hp coincide with the 

operator H k+2 M- Therefore, the eigenvalues A of problem (14) in the boson and fermion 
cases coincide with the corresponding eigenvalues of the operators H N (12). For the case 
in which the commutators between the operators b~(x, s) and b + (x, s), as well as the 
operators B~(x,y) and B + (x,y) are 1/N small, the asymptotics of the solutions of the 
problem (14), according to [15], are determined by the extremum points of the symbol 
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corresponding to problem (14). In the boson case, the symbol has the form 
H B [b*(.),b(-),B*(-),B(.)} 

/ eye ^ oo oc n 

= I E E ■ ■ ■ E / • • • / dxi ■ ■ ■ dXkdyi ■ ■ ■ dy2M x 

^ k,M=0 ' ' si=0 s k =0 J J 

xb*(x 1 ,s 1 ) b*(x k ,s k )B*(y 1 ,y 2 ) B*(y 2M _ 1 ,y 2M )H k+2M x 

x Symm (&(a;i,si) b(x k , s k )B(y 1 ,y 2 ) B(y 2M _ 1 , y 2M )) \ x 

xi...x k yi...y 2 M ^ ' J 

x I E J^Mn E • ■ • E / • • • / dx '^ ■ ■ ■ dx '^ ■ ■ ■ dy '™' x 

^ k',M'=0 ' ' s' n =0 s ' fe /=0 



x6*(xi,si) 6*(4„4)S*(yi,^) B*(y' 2MI _ 1 ,y , 2MI ) x (15) 

x Symm (&(ar'i,si) 4') 5 (^> 2/2) B (l/2M'-i,l/2M') • 

*i-*' t ,Vi-l^, V 7 J 

In the Fermi case, the symbol is expressed similarly, except that Symm in formula 
(15) is replaced by Asymm. 

The following identity for the symbol (15) holds in the Bose case: 

H B [b*(.)M-),B*(-),B(.)} = ^^, (16) 
where H, ps are the secondary quantized operators 



H = J dx$ + (x) ^ - -^—Aj$-(x) + 



2 

Here ps depends on the functions b(x, s), B(y,y'): 



+\ I I dxdyV(x,y)il) + (y)i> + (x)i> (y)ip (x). (17) 




OOOO s OO „ „ \ k 

p* = E E kWk + 2M)\ E / / dxdx ' b ^ °w> >< 

x ( // dyidy2B ( yi > y^ + (y^ + (y^ x 

x^JJ dy[dy' 2 B{y[,y' 2 )i-(y[)r(y' 2 )^ x 

x exp ( - / dz$ + (z)$~(z) j , (18) 



where ip + (x), ^ (x) are the Bose creation and annihilation operators, ordered according 
top Vick [26]. In the Fermi case we have a similar identity 

n F r(.),b(-),B*(-),B(.)]^-^0, 

where H, p F are the following secondary quantized operators: 



H — / / dxdyV(x,y)i/j + (x)i/j + (y)ip (y)ip (x) 




2 
and 

= EE fc!M!(A; + 2M)! ( / J dy 1 dy 2 B(y 1 ,y 2 m yi my 2 )j x 

oo OO „ „ 

■ ■ ■ / . . . / dxidx[ . . . dx k dx' k x 



S1=0 S fc =0 ' 



x6(xi, si)^*^!, si) b(x k , s k )b*(x k , s k ) x 

x^ + (a;i) ? + (x fc )Po?~(a4) x 

x( ffdy[dy 2 B(y[,y' 2 )r(y' 2 )r(y[)] ; (19) 



here, in the given case, ■?/' + (x), ^ (x) are the Fermi creation and annihilation operators 
and Pq is the projector on the vacuum vector of the fermionic Fock space. In general, for an 
arbitrary secondary quantized operator A, the symbol of the corresponding ultrasecondary 

quantized operator A is expressed [15] by the following formula: 

A B , F [b*(-),b(-),B*(-),B(.)} = S g P ^^ ) - 

In the space J 7 , we introduce [15] the ultrasecondary quantized operators for the number 
of particles: 

oo oo oo oo 



^ = EE< fc + 2M)ngr , N F = £ £ (Jfc + 2M)n^ mm - (20) 

fc=0 M=0 fc=0 M=0 

Correspondingly, in the boson case, the operator Nb has the form 

^= EE -^r E • • • E / • • • / • • • dx ^M x 

^ k=0 M=0 ' ' si=0 s fe =0 >/ 7 

x6*(xi,si) 6*(x fc ,Sfc)S*(x fc+ i,x fc+2 ) 5*(x fc+ 2M_i,a;fc+2Af) x 

x Symm (&(a;i,si) b(x k , s k )B(x k+1 , x k+2 ) ■ . . . 

Xl...Xk+2M ^ 

■B(x k+2 M-1, X k+2 M) ) f X 



OO OO -| oo oo 

x | E E E--- E / ••• d Zl ...d Zk , +2M , 

H'=0M'=0 ' ' s[=0 s' k ,=0 

•b*(zi, s[) b*(z k >, s' k ,)B*(z k / + i, z k i +2 ) B*(x k ' +2 M'-i, x k i +2 M>) x 

x Symm (b(z u s{) B(z k , +2M >-i, ^'+ 2 m')) [ • ( 21 ) 



Z l- Z fc'+2M' 

In the corresponding formula for the fermion case, Symm is replaced by Asymm. 



2 Symbol of the ultrasecondary quantized operator 



First of all let us note that the definition of the symbol given above does not reflect the the 
thermodynamical asymptotics, although this definition is in accord with the Bogolyubov- 
Dirak rule asserting that the creation and annihilation operators in the leading term of 
the asymptotics must be replaced by c-numbers. However, we shall say that the symbol 
obtained in this way is a pseudosymbol. Suppose that the operator H is of the form 

H = ^2 / • • • / d%\ • • • dxiip + (xi) . . . ip + (xi) x 
1=1 

i i 

(2 2 (9 \ ^ ^ 

x u . . .,xi\ -*7^r> • • • » )^~{ x i) ■ (22) 

Then, in the case of ultrasecondary quantization without the creation and annihilation 

operators of pairs of particles B ± (x, y) for the operators H and E defined above, we have 
the identity 

H = EA, (23) 
where A is the operator on the space J 7 , of the form 

L 00 



a = ^2 ^2 ' ' ' ^2 / • • • / dxi ■ ■ ■ dx ^ + i x ^ s • • • $ + ( x h so 



x 



1=1 si=0 s t =0 



1 1 



(2 2 d \ ^ ^ 

x u ...,xi\ •••■> ~ i g~rJ b ^ Xh Sl "> ■■■ b ~( x i)- ( 24 ) 

For the case in which the ultrasecondary quantization also takes into account the 
creation and annihilation operators for pairs of particles, the identity (23) remains valid, 

except that the operator A has a more complicated expression than the one in (24). For 
example, when we consider the operator H in the particular case of the formula (22) with 

L — 2, the operator A will have the form 



A = J2j dxb + (x,s){^- —Ajb~(x,s) + 



+ // dxdyB + (x,y)( - ^-{A x + A y ) )B-(x,y) + 



+ 2 Yl // dx dyV(x,y)b + (x, s 1 )b + (y, s 2 )b~(y, s 2 )b~(x, si) + 



S 1= 0s 2 =0 
OO 



s=0 



+ / / / dx dy dz(V(x,y) + V(x,z))xb + (x,s)B + (y,z)B~(y,z)b~(x,s) + 



+ I I dx dyV(x,y)B + (x,y)B (x, y) + - I I I I dx dy dz dwV (x, y)B + (x, y)B + (z, w) x 




2 

x^£-(y,w)£-(x,^ + £"(«;, 

If, in the expressions (24), (25), we replace the operators B ± (x,y),b ± (x,y) by c-numbers, 
we obtain the symbol corresponding to the asymptotics in the thermodynamical limit. 

in 



Consider the system of N identical bosons of mass m locate in the three-dimensional 
rectangle T with side lengths L±, L 2 , and L 2 . We assume that the bosons interact with 
interaction potential 

V(N^ 3 (x-y)), (26) 

where V(£) is a finite even function, x, y are the coordinates of the boson in the rectangle 
T. The boundary conditions along L x will be assumed periodic, while along L 2 we impose 
the condition of having zero derivatives. Note that the potential (26) depends on N in 
such a way that the radius of the interaction potential decreases with the increase of the 
number of particles N, but so that on the average the number of particles which interact 
with a fixed particle remains constant. 

Under ultrasecondary quantization for pairs of bosons in the considered system, we 

obtain the (ultrasecondary quantized) operator H whose explicit form was presented 
above. As was explained previously, for this ultrasecondary quantized operator we have 
the identity 

H = EA, (27) 

where E is the ultrasecondary quantized unit operator, while A is an operator in the space 
of ultrasecondary quantization. It is easy to verify that the following operator satisfies an 
identity of the form (27): 

A = J J dxdyB + (x,y) ^-^-(A x + A y ) + V (N 1/3 (x - y))^j B~(x,y) + 

+2 j dxdydx'dy'V(N l l\x-y))B + {x,y)B + {x\y')B~{x,x')B-(y,y'), (28) 

where B + (x,y) and B~(x,y) are respectively the boson creation and annihilation op- 
erators for pairs of particles in the Fock space of ultrasecondary quantization. By the 
identity (27), in order to find the asymptotics of the spectrum of the boson system under 
consideration in the limit as N — > oo we must find the corresponding asymptotics for the 
operator (28). 

Since the function (26) of multiplication by N in the limit as N — > oo converges in 
the weak sense to the Dirak delta function, it follows that the second summand in the 
operator (28) in this limiting case contains the small parameter 1/N. This means that, 
in order to find the asymptotics of the eigenvalues and the eigenf unctions of the operator 
A one can apply the quasiclassical methods developed in [13]. The asymptotics of the 
eigenvalues and the eigenf unctions is determined by the symbol of the operator (28); 
this symbol is called the true symbol of the ultrasecondary quantized problem. The true 
symbol corresponding to the operator (28) is the following functional, defined for any 
pair of functions § + (x, y) and y): 

H [$+(•),$(•)] = jj dxdy$ + (x,y) (-^(A* + A„)) y) + 

+2 J dxdydx'dy' (NV (N 1/3 (x - y))) $ + (x, y)$ + (x', y')<S>(x, x')$(y, y'). (29) 

From the invariance of the number of particles in the system with the functions Q + (x, y) 
and y), we obtain the condition 
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J I dxdy$ + (x, y)Q{x, y) = i (30) 

According to the asymptotic methods (see [13]), to each solution of the system of 
equations 

"•<*•»> = «3b- m+ ^ = Jh) (31) 

which also satisfies condition (30), there corresponds, in the limit as iV — > oo, the asymp- 
totic series of eigenvalues and eigenfunctions of the operator (28). From the expression 
for the true symbol (29), it follows that the system of equations (31) can be written in 
the form 

h 2 

m(x, y) = - — (A x + Ay)$(x, y)+ 

+ JJ dx'dy' (NV (N 1/3 (x - y)) + NV (N 1/3 (x' - y'))) $ + (rr', y')$(x, x')$(y, y'), 

h 2 

^+(x, y) = - — (A x + Ay)$ + (x, y) + 

+ JJ dx'dy' (NV (iV 1/3 ( x — x')) + NV (A^ 1 / 3 (y — y'))) y')& + (x, x')& + (y, y'). 

Let v q be the Fourier coefficients on the box (L 1; L 2 , L 2 ) of the potential NV(N 1 ^ 3 x): 

v q = — L J e~ iqx NV (^N(x^j dx. (32) 
The exact solution of system (32) is given by the following functions: 



L_ e -ik 1 ( X +y) C0S ^ x _ y ^. (33) 



2 



$fc 1; fc 2 



l-J^tPk^^e^-rt (34) 



with eigenvalues 



h 2 

n = — ( k\ + k\) + v + v 2k2 , (35) 
m 



where the function (p k2 j is given by 

bi , 1 



/i 2 /m(/ 2 - £|) - ( Uo + w 2fc2 ) 



^fc 2 ,fc 2 



W;_ fc2 + Vi +k2 2 

The pair of vectors k\, k 2 plays the role of parameters indexing the various solutions 
of this system. The vector hki/m expresses the velocity of the flow of the boson system 
along the capillary. The vector k 2 is the wave vector of the transversal mode. 
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Note that as |/| — > oo we have bi — > oo, since 



N = 7^72 / e^V(Ode < / ^(OK < 7TI / < vo, (36) 

-^1-^2 •/ ATT lj \ lj 2 J NT ij \ lj 2 -R 3 

so that the series (34) converges absolutely. 

Let us split the series (34) into two parts for / < N 1 ^ and for I > iV 1//6 . The first sum 
converges as N — > oo with precision up to (iV -1 / 6 ) to the value 



2mVb 

This follows easily from the substitution = £ in equation (32). 

The second part of the sum tends to zero by (??) up to 0{N^ 1 ^). Therefore in the 
limit as N — > oo, the system of equations (32) under the additional condition (30) 
possesses, for k\ = 0, the following family of solutions: 



®t y) = jji cos ( k (x - v)) 



^k{x,y) = -^— 1 S2^ k) ie^(il{x-y)), (38) 



where k and I are three-dimensional vectors of the form 

n 2 n 3 * 



27T ( 0, , , 

here n 2 , n 3 are integers and the terms ipk,i in formula (27) have the form 



= 2^ ^ " + V » ± V (£< P " P > + V ») " V ° 2 I • < 39 » 

where the plus sign is chosen when I 2 > k 2 , the minus sign, when I 2 < k 2 , and, finally, Vo 
here, as before, stands for the quantity 



V " = TJ 2 J dxV{x) > (40) 

in which the integral is taken over the space R 3 . 

The vector k in formulas (38) plays the role of a parameter indexing the various 
solutions of equations (32), (30). The solutions (38) are standing waves and correspond 
to series for which there is no flow. 

The leading term of the asymptotics of the eigenvalues of the series corresponding 
to the solution of (33) (34) is equal to the value of the symbol (29) on these functions 
multiplied by iV: 
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The asymptotics of the eigenvalues and eigenfunctions, in particular the terms that 
follow Ek lt k 2 , are determined not only by the system (32), but also by the solution of the 
variational system for the Hamiltonian system of equations. The system of variational 
equations for (32) has the form: 

(Q - A) F(x, y) = ~(A X + A y )F(x, y)+ 
2m 

+ 2N JJ dx'dy' (v(^N(x-y)) +V(^N(x'-y'))^j (G(x' ,y')<$>(x,x')<$>(y,y') + 
+ $ + (x', y')F(x, x')$(y, y') + <$> + {x', y')$(x, x')F(y, y')) , 
(n + A) G(x, y) = -^(A x + A y )G(x, y)+ 

+ 2N JJ dx'dy' (y(^N(x-x')) + V(^N(y-y'))^ (F(x',y')$ + (x,x')$ + (y,y') + 
+ <f>(x', y')G(x, x')<f> + (y, y') + $(x', y')<$> + (x, x')G(y, y')) . 

(42) 

In order to find the spectrum of the quasiparticles, one must distinguish, among all the 
solutions of the variational system, those which satisfy the selection rule 

JJdxdy (G*(x,y)G(x,y) - F*(x,y)F(x,y)) > 0. 

If k 2 — 0, then the asymptotic series corresponding to this solution is the Bogolyubov 
series with flow velocity hk\/m. The spectrum of quasiparticles of this series is expressed 
by the well known formula 



Consider the case k 2 ^ 0. Substituting the solution (33), (34) into (42) and tak- 
ing symmetry into account, we obtain the following expression for the solution of the 
variational system: 

Gi(x, y) = ^ ^ exp(i(A;i + k 2 )x + i(ki + l)y) + exp(i(A;i + k 2 )y + i(h + l)x) j + 
- ( exp(i(/ci - k 2 )x + i(ki + 2k 2 + l)y) + 



+ exp(i(A; 1 - k 2 )y + i(k t + 2k 2 + l)xf) , (44) 

Fi(x,y) = ^-(^exp(i(ki + k 2 )x + i(k\ + l)y) + exp(i(A;i + k 2 )y + i(k\ + /)x)j + 

+^ ( exp(i(A; 1 - fc 2 )x + i(h + 2k 2 + l)y) + 

+ exp(i(fci - k 2 )y + i(h + 2k 2 + l)xfj + 

+ Yl ^[exp(i(k 1 + k 2 + l-l')x + i(k 1 + l')y) + 

l'^l,l+2k 2 

+ exp (i(fci + A; 2 + / - + + Z')a:) ) , 
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where I ^ —k?, while the numerical coefficients uij, U2,i, vij, V2,i, wiy are determined 
by an infinite system of equations. This system contains a closed subsystem consisting 
of four equations for the coefficients Uij, 112,1, v\,i, V2,i, which can be written in standard 
form 

AX = MX, 



(45) 



where 



here X is the column vector 



A = A 



X 



h 2 



fci(fc 2 + Z); 



M is the matrix 



M = 



Bi 



\ V ((Pk 2 ,l + Vk 2 ,l+2k 2 ) 

and the B t have the form 

K 2 



m 

( \ 

u 2 ,i 

V l,l 



V 

Bl + 2k 2 

Vo i.Vk 2 ,l + Vk 2 ,l+2k 2 ) ~B, 
,l+2k 2 



-V 





\ 





-Vo 




-B t 


-Vo 




-V 


— Bl + 2k 2 


/ 



A 9 



V (pk 2 



The matrix M can be written in block form 

M = 



C -V E 
D -C 



C 



D 



2V (fik 2 ,i 

Vq (<fk 2 ,l + fk 2 ,l+2k 2 j 



Vo {<Pk 2 ,l + <fk 2 ,l+2k 2 ) 
,l+2k 2 



constituted by the following 2x2 matrices: 

Bi Vo 
Vo B t+ 2k 2 

and by the unit matrix E. 

This expression allows us to rewrite equation (45) 

\ 2 u = {C 2 - V D)u 

for the column vector u — \ Ul ' 1 

The selection rule (44) for the variational system of equations, which has the form 

+ u h u ^ ~ v u v u - v h V2 i > °- 

specifies the following eigenvalues: 

h 2 

><iMM,i = k i( k 2 + l)± 



h 2 

2m 



I* 
2 



I 1 
_i 

2 



H + - + - J4C 2 + >?) ^(i 2 + >i-^)\/(> 2 ->i) + 



AmVo 



h 2 



(46) 



^2,fci,*a,J = + l)± 

m 



h 2 



± + l l + I - k 2 (P + l\) - ±(P + 2k 2 ) ] J(P - l\) + (^), 
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where l\ — I + 2k 2 and / ^ k 2 . In the formula for Ai^^, we choose the plus sign if 
I 2 > k 2 , and the minus sign if I 2 < k 2 , whereas in the formula for \ 2 ,ki,k2,i, we take the 
plus sign if l\ > k 2 , and minus, if l\ < k\. 

Formulas (46) for I ^ —k 2 determine the spectrum of the quasiparticles of the series 
corresponding to the solutions of (33), (34). From the explicit form (46) it follows that 
there are negative values in the spectrum of the quasiparticles. Therefore, the series 
corresponding to the solution (33), (34) for k 2 ^ is not metastable. We shall assume 
in what follows that L\ ^> L 2 . Consider the Bogolyubov series corresponding to the flow 
along the capillary with velocity hko/m, where k = 2ir(ni/Li, 0, 0). For the system of 
bosons under consideration the leading term of the asymptotics of the eigenvalues equals 

N (47) 

Now let us assume that the relationship between L\ and L 2 is such that there exists 
a pair of vectors k±, k = 2n(0,n 2 /L 2 ,n3/L 2 ) for which the corresponding value (41) is 
exactly equal to (47). This means that there may be resonance between the superfluid 
states of the Bogolyubov series and the nonflowing states of the metastable series corre- 
sponding to (33), (34). If Li is very large, resonance is also possible for the case in which 
the value (41) is close to (47) but does not necessarily coincide with it. 

The existence of such a resonance makes possible the passage from the superfluid state 
to the metastable one, from which the system drops to the lower energy level, which means 
that superfluidity is lost. The minimal energy arising in the metastable series corresponds 
to the case k = 27r(0, l/L 2 , 0) and, according to formula (41), equals 

Comparing (48) with (47), we see that resonance is impossible if the absolute value of 
the flow velocity v is less than a certain number 

. . 2nh h , 4 . 

h s szr (49) 

The quantity v c (L 2 ) in the right-hand side of inequality (49) increases as the cross 
section of the capillary L 2 decreases. If L 2 is less than 27rh/ (mv C L), where v c l is the critical 
Landau velocity determined from (43), then v c (L 2 ) is greater than the critical Landau 
velocity and, accordingly, in this case superfluidity disappears when the Landau velocity is 
reached. However, when L 2 increases, the resonance between the flowing and metastable 
states arises before the Landau velocity is reached, and this explains the dependance of 
the critical velocity on the size of the cross section of the capillary. 

Since the Landau curve and the Bogolyubov quasiparticles, as we showed in [10], do 
not change in the classical limit, it follows that, under the condition v < h/(mL 2 ), where 
L 2 is the diameter of the nanotube, we have the following mathematical fact: the classical 
liquid in the nanotube must be superfluid. 

Thus we have obtained a series of quasiparticles that differ from the Bogolyubov series. 
This series corresponds to a "standing wave" across the tube, while the Bogolyubov wave 
is a "running wave" along the tube. 

The question of the dependence of these series on the temperature arises. As was shown 
by the author previously, the Bogolyubov series depends on the temperature, which leads 
to the dependence of the Landau criterion (critical level) and to a new phase transfer, 
which I called phase transfer of the zeroth kind. 
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It turns out that for the "transverse series" presented here and in [21, 22], the leading 
term of the asymptotics as iV — > oo does not depend on the temperature. This follows from 
the fact the operator A in the corresponding (L, A)-pair which determines the dependence 
on temperature in its leading term "vanishes" as iV — > oo in the sense that A 2 = 0. 

Thus, for a sufficiently narrow tube, the phase transfer is determined by the Bo- 
golyubov series only, so that here we have a phase transfer of the zeroth kind. 

Remark 

In our problem, the asymptotic decomposition is considered with respect to two pa- 
rameters: N — > oo and h — > (the quasiclassical limit). How are they related? In 
the 1947 Bogolyubov paper, there is no external potential, and so the problem is solved 
exactly. This means that, just as for linear equations with constant coefficients, the qua- 
siclassical solution coincides with the exact one. But what if we introduce a potential? 
In that case, as we mentioned in the introduction, the quasiclassical approximation yields 
quasiparticles that preserve superfluidity (see [8] formula (24) and [9]). 

However, first of all, the following conditions for the external potential must hold: in 
the classical problem of a self-consistent field, there must exist an invariant Lagrangian 
manifold up to "atomic precision", i.e., the potential must be "absolutely precise" and no 
noise of order h is acceptable. 

If this is not the case, the Bogolyubov spectrum in the three-dimensional case is 
destroyed (see [6] and [15]) and the superfluidity effect disappears. 

Secondly, the potential must change very slowly, so that its change on a finite interval 
be less than a value of order 1/N. Such a "nanoprecision" and the presence of such drastic 
conditions is possible only if nanotechnologies are used. 

I surmise that a verifying experiment should be performed by using a neutral gas such 
as argon. 

Let us note in conclusion that high temperature superfluidity, which was theoretically 
established by the author in 2004 [27, 28, 29, 30, 31], is essentially the superfluidity of 
liquids in nanotubes. 

I would like to express my gratitude to D.S.Golikov for recalculating and verifying all 
the formulas. 
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